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Abstract 



T— I jThis paper studies composite quantum systems, like atom-cavity systems and coupled optical resonators, in the absence of external 
driving by resorting to methods from quantum field theory. Going beyond the rotating wave approximation, it is shown that 
the usually neglected counter-rotating part of the Hamiltonian relates to the entropy operator and generates an irreversible time 



evolution. The vacuum state of the system is shown to evolve into a generalized coherent state exhibiting entanglement of the modes 
in which the counter-rotating terms are expressed. Possible consequences at observational level in quantum optics experiments are 
currently under study. 
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I 1. Introduction 

This article concerns one of the most useful approximations 

ajn quantum optics and atomic physics, namely the so-called 
I I rotating wave approximation (RWA) 1^1. Up to now, there is 
jn general very good agreement between experimental findings 
and theoretical predictions based on this approximation. Differ- 
>^ ent reasons for the validity of the RWA are given in the litera- 
ture. Most authors argue with time scale separation. Indeed, in 
jhe presence of sufficiently weak resonant interactions (Uke res- 
onant laser driving of optical transitions) it is possible to move 
' into an interaction picture, where the counter-rotating terms os- 
cillate very rapidly. Their contribution to the time evolution of 
the system hence remains negligible when compared with the 
effect of the non-rotating terms |2]. Other authors apply the 
RWA in order to preserve quantum numbers and energy yj, 0] 
and the validity of the so-called two-level approximation ISO. 
, However, recent trends in experimental quantum optics and 
atomic physics aim at the realisation of miniaturised devices 
10, 0] with coupling constants which are many orders of mag- 
nitude larger than comparable, more classical designs |8]. As 
a result, the separation of the relevant time scales in a system 
might be reduced by many order of magnitude such that conse- 
quences of the normally neglected counter-rotating terms in the 
system Hamiltonian become observable. Systematic studies of 
the qualitatively different parameter regime, i.e. beyond the ro- 
tating wave approximation, are currently becoming feasible in 
quantum optics experiments. An example, where the counter- 
rotating terms are already routinely taken into account, is the 
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calculation of temperature limits in laser and in cavity cooling 
experiments jgl lToll . 

But there are other situations in which the RWA should not 
be applied. Also in certain situations far away from resonance, 
the counter-rotating terms should not be neglected ifTHl . An 
example is discussed in Ref. by Hegerfeldt, who showed 
that the interaction between two atoms and the free radiation 
field, when treated exactly, can result in a small violation of 
Einstein's causality. Zheng et al. |13] also avoided the RWA 
and predicted corrections to the spontaneous decay rate of a 
single atom at very short times. Recently, Werlang et al. 
and the current authors f\3\ pointed out that it might be possi- 
ble to obtain photons by simply placing atoms inside an optical 
cavity. A similar energy concentrating effect might contribute 
significantly to the observed very high temperatures in sonolu- 
minescence experiments [161. 

What Refs. US [H Hlil [3 [Hi have in common is 
that all of them consider open quantum systems, i.e. systems 
coupled with environment or with other systems (acting as en- 
vironment). In such open systems energy is no longer necessar- 
ily a preserved quantity and time evolution is irreversible. For 
example, Refs. |14, 15] predict an energy concentrating mech- 
anism in coupled atom-cavity systems, even in the absence of 
external driving. This might seem unphysical but does not con- 
stitute a violation of the laws of thermodynamics, as long as 
the predicted changes in the free energy are accompanied by 
respective changes in entropy. 

In this paper, we consider a composite quantum system 
which is much simpler but nevertheless closely related to the 
one considered in Refs. ifTl. [Tsll . Resorting to quantum field 
theory (QFT) methods, we show that the two linearly-coupled 
bosonic reservoirs shown in Fig. [1] are characterised by a non- 
trivial entropy operator caused by the coupling of field modes 
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Figure 1: Schematic views of two possible experimental realisations of two 
linearly-coupled bosonic reservoirs, (a) The experimental setup considered here 
could consist of a large number of strongly-confined two-level atoms inside the 
free radiation field or inside a relatively large optical resonator with a finite 
bandwidth, (b) The same Hamiltonian is obtained for two coupled optical res- 
onators with overlapping field modes. 

and generating an irreversible time evolution. We show that 
the origin of the dissipative dynamics are indeed the usually 
neglected counter-rotating terms in the interaction between dif- 
ferent components of the system. 

In particular, we study the meaning and the effects of the 
RWA from the perspective of the space of the states of the sys- 
tem under study and show that when dealing with quantum 
fields, namely with infinitely many degrees of freedom, as it 
is necessary in the study of open systems, the counter-rotating 
terms cannot be neglected. This can be understood in general 
terms by observing that, as far as one limits himself to systems 
with finite numbers of degrees of freedom, the von Neumann 
theorem in quantum mechanics (QM) guarantees that the rep- 
resentations of the canonical (anti-) commutation rules (CCR) 
(the Hilbert spaces of the states) are unitarily equivalent and 
therefore physically equivalent [ ITillllllSl- Thus, in QM there 
is no room for non-unitary transformations. One does not need 
to wonder about the choice of the representation in which the 
system dynamics, i.e. the Lagrangian or the Hamiltonian opera- 
tor, is realized since all of them have the same physical content; 
they are, indeed, equivalent up to a unitary transformation. 

A different situation occurs when one deals with infinitely 
many degrees of freedom, as it happens when field operators 
are considered, i.e. in QFT. In such a case, the von Neumann 
theorem does not hold anymore and infinitely many unitarily 
inequivalent representations of the CCR exist |17, 18, 19]. In 
this case, the choice which one among these representations 
to adopt, i.e. where to realize the system dynamics, might be 
of crucial physical relevance. For example, one might real- 
ize the system dynamics in a representation where the symme- 
try properties of the ground state are the same as those of the 
field equations, or one might instead use a representation where 
spontaneous symmetry breakdown occurs llisl I20I1 . Also, one 
might choose a representation where the ground state is pre- 
served under time evolution, or, instead, a representation where 
it changes (as in unstable systems) under time translation trans- 
formations |21]; and so on. 

The formal and physical significance of the unitarily inequiv- 
alence among representations is that the vacuum state in each 
of them cannot be expressed in terms of the vacua of other rep- 
resentations. Thus, for example, the vacuum of a metal in the 
superconductive phase cannot be expressed in terms of the vac- 
uum of the (same) metal in the "normal" phase. In phenomena 



related with unitarily inequivalence, a dominant role is typically 
played by "weak" couplings (e.g. in the field theories, the 
order parameter, say a, which specifies the vacuum, is given 
by oc \irp-\/A, where A is the coupling constant and is 
the (negative) squared mass). These weak coupling effects can- 
not be studied in a perturbative expansion around the vanishing 
value of the coupling constant (e.g. in the A<p'^ model, the or- 
der parameter a cannot be defined at A - 0). As we will show, 
deciding whether or not to apply the RWA can constitute a sim- 
ilar delicate type of problem which, as said, cannot always be 
ignored when dealing with quantum fields. 

Our calculations in the interaction picture reveal the under- 
lying degrees of freedom which are involved in the generation 
of entropy and non-unitary time evolution in the setup shown 
in Fig. [1] We show that the initial vacuum state of the system, 

1. e. the state with no population in either mode, evolves in time 
into a generalized SU( 1 , 1 ) coherent state which is orthogonal to 
the initial state and exhibits entanglement of the modes in which 
the counter-rotating terms are expressed. The representation of 
the system then requires a different Hilbert space at any mo- 
ment in time. This constant changing from one state space into 
another makes the dynamics of the system irreversible. When 
applying the RWA, the time evolution of the system remains 
unitary and the vacuum state of the system is always the same. 

There are five sections in this paper The theoretical model is 
introduced in Section|2] In Section|3]we show that the counter- 
rotating terms in the Hamiltonian of this system are related to 
the entropy operator and a detailed proof of ureversible time 
evolution is presented. In Section |4] we study the free energy 
and the entanglement of the vacuum state. We show that the 
free energy is minimized at each time t of the evolution of the 
vacuum state. Our discussion thus leads us to conclude that in 
conditions far off resonance the counter-rotating terms in the 
Hamiltonian are related to entropy, free energy and entangle- 
ment. These results are summarized in Section |5] In the Ap- 
pendix we report some mathematical formulas used in the text. 

2. Theoretical model and time-dependence of the vacuum 
state 

For the sake of simplicity and concreteness, it is convenient 
to refer to a specific model (which is widely used in quantum 
optics). Our discussion, however, can be extended to other 
models. We thus consider for example the familiar Hamilto- 
nian H for an ensemble of tightly-confined two-level systems 
(e.g. atoms)coupled to a radiation field (cf. Fig. [Ha)). It can be 
written as llj] : 

H^Hq + , (1) 
where 

k / ^ 

^^int ^'^YjYj ("k ~ "k) (o-t + o-r) . (2) 

k / 

The cr,'s are the two-level system SU(2) spin-like operators, 
flk and are the boson annihilation and creation operators of 



photon modes k. The corresponding characteristic frequencies 
are a>o and cjk- In the dipole approximation, the (real) atom- 
field coupling constants in S is given by 



8k 



IheoV 



1/2 



doi ■ e 



(3) 



where doi is the dipole vector, e is the polarization vector, and 
V is the volume. For simplicity we assume that the atoms are 
well localised within an optical domain and e'*""' = 1 for all 
wave vectors k and atomic positions r,. Moreover, all atoms 
have the same dipole moment doi. This is why we consider 
only one photon polarisation and why the do not depend on 

Suppose |0)j denotes the vacuum for the two-level atoms, 
|0)r the vacuum state for the radiation field, and |0) = |0.s, Or) = 
|0), ® |0),.. Then cr-|0) = = flklO). 

Alternatively, one may consider a system of radiation field 
modes at and interacting with another set of reservoir (or 
cavity) field modes with boson annihilation and creation oper- 
ators and b\^, as illustrated in Fig. [Ub)- The Hamiltonian H 
can again be written as in ([TJ. Adopting a notation analogous 
to the one above, we now have b\f\Qi) = = flklO). The Hamil- 
tonians Hq and //int are then given by 



k k 



(4) 



where the dipole approximation has been used again and where 
the Qjt are the frequencies of the b]^ modes. Note that the (bo- 
son) operators at and fet commute ([ak, ^] - 0, and all other 
commutators are zero). The commutation relations for them are 



[flk, Cq] - ^kq , [^k, ^q] - ^kq ■ 



(5) 



All other commutators vanish. 

It should be recalled here that the boson operators b]^ and b^ 
can be related to the spin-like operators cr"^ by considering the 
ensemble of two-level systems (atoms) for large A^. This is 
discussed in detail in Ref. ll22tl . Since this derivation is outside 
the task of the present paper we omit to repeat it here again. 
We only recall that in the large limit (A^ — > oo) the Weyl- 
Heisenberg algebra is obtained as the contraction of the su(2) 
algebra for the cr's |22]. As a result, in the large A^-limit, both 
situations illustrated in Fig. [1] can be described by Hq and Hint 
given in Eq. For definiteness, we consider these in the fol- 
lowing. 

Since we are here especially interested in the effect of the 
counter-rotating terms in the Hamiltonian Hint, it is convenient 
to write 



^^int - Hjc + HcR 



(6) 



with 



k 



(7) 



The Hamiltonian Hjc is the usual Jaynes-Cummings Hamilto- 
nian ||23ll and the part that usually survives in the RWA, while 
HcR is the counter-rotating term part. There exists an interac- 
tion picture in which these counter-rotating terms are fast os- 
cillating terms and they are "therefore" neglected. The RWA 
consists in fact in neglecting those terms whose (antiresonant) 
frequencies (^ (L)k,+ = oJk + ^k) are far off from the resonance 
condition {clJu- = ojit - Q.^ ~ 0). 

From Eqs. © and O we see that the vacuum |0) is annihi- 
lated by //() and Hjc but not by Hcr'- 



Hq |0) = 0, H,c 10) = 0, Hcr |0> ^ . 



(8) 



This means that the vacuum of the theory is not invariant under 
time-translation unless the RWA is applied. This is an inter- 
esting feature from the mathematical point of view as well as 
from the physical point of view. It reminds us of the mecha- 
nism of spontaneous breakdown of symmetry. In the present 
case, the spontaneously broken symmetry is the one of time- 
translation: the dynamics is invariant under time-translation 
since obviously [H, H] = 0; however, the vacuum is not in- 
variant since e^"^^^ |0) |0) due to Eq. ([S). So, let us consider 
this feature in more detail. 

For notational simplicity, we focus our attention on only one 
of the k modes and omit the suffix k in the following. At the 
end we will recover them. In the Appendix we introduce the 
generators 7,- and /,■ (/ - 1, 2, 3) of the SU(1,1) and SU(2) group, 
respectively. In terms of these generators, the Hamiltonian H - 
Hq + Hjc + Hcr with Hq, Hjc and Hcr given by (|4]i and (|7]l can 
be written as 



■2gl2- 2gJ2 



for each single k mode. For simplicity we have set h 
Notice that we have 

Hcr -272 , 



(9) 



(10) 



i.e. Hcr is, apart from the coupling factor, nothing but the J2 
generator of the SU(1,1) group (cf. Eq. (IA.2I )). 

In order to study the effects of Hcr (cf. Eq. ([8])) on the vac- 
uum state, one could directly compute e""^cR |0). However, it 
is instructive to see how the time evolution operator e^"" it- 
self, with H given by Eq. acts on the vacuum. To do that 
we rotate the state e^"^ |0) into a frame which is more conve- 
nient for our study. This requires two successive rotations of 
e^"" |0) which are induced by the generators /i and K2, respec- 



tively, with K2 being the squeezing generator 112411 given in the 
Appendix (cf. Eq. ( IA.17l i). We thus consider the following ro- 
tations 



-irH 



-itH(e) 



10) 
10) 



10) = e-''«(»> 10) , 



The rotation angles 6 and a are fixed by the theory parameters 
as 



2^ 



tan 6 - — , tanh a 



a)+A 



(12) 



respectively, with A and B given in the Appendix 
(cf. Eq. iAl2\) ). H(e) and H(e,a) are given by 
H(0) = e""'H e-""' and H(e,a) = e'"'^'-H{0) e-''"^\ re- 
spectively (cf. Eqs. (IA.25I ) and ( IA.26I )). We also have 
e'Wi |0> =_I0). The state lOCo-)) = e'"*^- 10) is the squeezed 
vacuum [24]. 

In the corresponding interaction representation, the time evo- 
lution of the vacuum is solely controlled by the interaction 
Hamiltonian in the interaction representation. There the a and h 
operators and their Hermitian conjugates carry their respective 
time dependence, a{t) and h{t) etc. Thus in the interaction rep- 
resentation with respect to the free Hamiltonian (cf. Eq. ( IA.26b ) 



(13) 



with £ = BjlA, the interaction Hamiltonian Hy^^{6, a) is (up to 
a constant term) given by 

i^:(.,.) = ^/3-^/..44^,. (14) 



The operator Ki is given in the Appendix (cf. Eq. ( IA.16I )). The 
operators h, J2, and K\ have to be understood to be in the in- 
teraction representation (i.e. they are made out of a(t) and b(t) 
etc.). 

We now observe that 



P<0(a)|, 



-itHZXe,a) 



10(a))" 



(15) 



{Q\e-''"^2 e-'^^'^r.^^''"^ e-'^'lO) 



= {<d\e-'""»''^\0) , 

where 

ffne,a)^Ah{t)-2YJ2{t) 



(16) 



with F = g a>-/A. Here we explicitly write /3(f) and 72(0 to 
remind the reader that we are using the interaction representa- 
tioiiii. Since 72(0 commutes with 73(0 (cf. Eq. jA.lll i). and 
73(0 |0) = 0, Eqs. O and (QSll yield 



'P<0(a)|e-'«'''«>|0(a))'P = (0|e''2"^«|0) . 
By using 

^m,h-^)^it2rj2(t)^-it6(j,-k) ^ ^it2Tj2 

and 

^±m,h-h^Qy = |0> , 

we finally obtain 

(0|e''2"2«|0) = <0|e''2"^|0) = (0|0(0) ■ (20) 

Here the notation |0(0> = e"^^^-\Q) has been used. Note that 
F ^ in the resonant condition limit w_ — > 0. From Eq. ( fTSI ) 
we see that the RWA is automatically impUed in such a Umit. 



(17) 



(18) 



(19) 



' It is interesting that the Hamiltonian of the quantum damped oscillator has 
the same form as Eq. )16t . See Ref. (2lh . 



Eq. ( l20b shows that the time evolution of the vacuum is in gen- 
eral non-trivial and indeed generated by an operator propor- 
tional to 72. 

In the next Section, we calculate an explicit expression for 
(010(0) by resorting to the well-known results of QFT iflslEl. 
l25n . Our discussion will show that the counter-rotating part 
of the system Hamiltonian, i.e. 7/cr, is related to irreversible 
time evolution, entropy, entanglement, and free energy. To our 
knowledge, these properties of Hcs. have so far been ignored. 
In conditions far off the resonance or in the absence of external 



driving they may produce observational consequences 11151 \\m . 



3. //cR, irreversible time evolution and entropy 

In order to compute the quantity (0|0(0) explicitly, we start 
by denoting the set of simultaneous eigenvectors of a ' a and b'^b 
by Wua, nh)]. The corresponding eigenvalues and nt are non- 
negative integers M25[l . Expressing 73 and (73 - 5), which in 
SU(1,1) form a complete set of commuting operators, in this 
basis, we find 

h\h m) = M m) , j = ^(Wfl - tih) , 

■^3-^1 1;, m) = m\j, m) , m = ^(n^ + iib) . (21) 



Here m > \i\ > since ria and «/, are non-negative 125|]. We re- 
mark that once the eigenvalue of 73 is set to be definite positive 
by boundary condition, then it remains constant under the time 
evolution induced by 72 (i.e. by 77cr), as it must be since 72 and 
/3 commute (73 is indeed the SU(1,1) Casimir operator). Also 
note that the original vacuum |0) of the system is actually the 
eigenstate of 73 and (73-5) associated with the zero eigenvalues 
7 = and m = 0, i.e. |0) = \n„ ^ 0, n,, ^ 0). 

We are now ready to compute the explicit expression for 
|0(0)- By use of the "normal form" of the operator e"^^^'^ (see 
e.g. chapter 4 of ref. |25] and refs. ifisl u\\\ ) we obtain the 



vacuum at time f 



|0(0) = e''^"-|0) 



1 



cosh(rO 

At each time t, |0(0) is a normalized state, 
(0(010(0) = 1 , Vf . 



exp(tanh(F07+)|0) . (22) 



(23) 



In the limit t 00, the vacuum |0(0) becomes orthogonal to 
the original vacuum |0) = |0(f - 0)). Indeed we obtain: 



(0|0(0) = exp (- In cosh FO ^0 . 



(24) 



The vacuum instability shown in Eq. (l24l) has to be expected on 
the basis of physical intuition since Hcr, being related with fast 
oscillating terms, introduces transient phenomena. These are of 
dissipative nature, their time evolution being controlled by e^^' 
for large t, as Eq. (l24l) shows. 

The time evolution induced by the counter-rotating terms in 
the Hamiltonian is thus shown to be only well defined for fi- 
nite short time-intervals (f < 1/F). As f — » 00, the time evo- 
lution manifests itself as a non-unitary transformation which 



4 



leads out of the original Hilbert space whose vacuum is |0). 
This is clearly a pathology, since due to the von Neumann theo- 
rem there is no room in QM for the non-unitary time evolution 
expressed by Eq. (l24l i. However, far off resonance (i.e. for non- 
vanishing one cannot neglect the counter-rotating terms. In 
this case, it becomes unavoidable to study the system dynamics 
without performing the RWA. 

We are thus led to explore the possibility to formulate our 
problem in the framework of QFT [21], where infinitely many 
unitarily inequivalent representations exist. In order to do that, 
we first restore the suffices k and write |0(f)) as 



10(0) = n 



1 



cosh (Fi?) 



exp(tanh(rif)A,k) |0> 



(25) 



where s Uk-jAk- Eq. (I25l l is a formal relation holding 
for finite volume V. As customary in QFT, one works at fi- 
nite volume and the limit V ^ oo is taken only at the end of 
the computation. The state |0(f)) is again a normalized state, 
(0(0|0(0) = 1, at each time t. In fact, it is an SU(1,1) general- 
ized coherent state f23\, i.e. a two mode Glauber-type coherent 
state. Eq. is now replaced by 



(010(0) = exp 



- ^ In cosh (F^O 



(26) 



which again exhibits non-unitary irreversible time evolution. 

In QFT, however, we have to consider infinitely many de- 
grees of freedom. Thus, by using the continuous limit relation 
^^ W(27r)^ / cf'k, we obtain 



<0|0(0) ^0 Vf, 
(O(f')IO(O) V f , f' with t + t' , 



(27) 



provided that Jc/^^A: lncosh(F/.0 is finite and positive. The 
meaning of Eq. (l27l i is that the representation at a given time t is 
unitarily inequivalent to the representation at any different time 
t' tm the infinite volume limit: The system spans a whole set 
of unitarily inequivalent representations as time evolves. Each 
of them is labeled by diff'erent values of t. The occurrence of 
such a phenomenon is possible in QFT where infinitely many 
unitarily inequivalent representations exist. Time evolution is 
thus described in terms of "phase transitions" among the repre- 
sentations, or "trajectories" in the space of the representations. 

The above calculations show that the generator of such a non- 
unitary time evolution is the counter-rotating term proportional 
to J2 = 2k ■^2,k- We now show that J2 is associated with the 
entropy operator 

The vacuum state |0(0) given by Eq. ( |25] | can also be written 
as filial 



10(0) = exp j -^s\ J) = exp l-^Sb )| J) . 



Here 



I J) = exp 



V k 



10) 



(28) 



(29) 



is a not normalizable vector |[ii|2ll]and5„ is given by 

Sa = - \\l a^Ok In sinh^^Fi?) - 0^0^ In cosh^^Fi?) > . (30) 

k 

Sb has the same expression with and b^^ replacing and a^, 
respectively. In the following, we write S for either Sa or Sb. 
It is not difficult to recognize that S is the entropy ifisll . Indeed, 
the state |0(0) can be written in the form 



\0(t))^j]^|w'„\n,n} 



«=o 

where n denotes the set {nt ) and 
sinh2"'(Fif) 
^ cosh2("'+i>(F^f) 
with 



< W„ < 1 and ^ W„ = 1 . 

n=0 

We have 

+00 

{0{t)\S\0(t)} ^ -Y^Wn(t)lnW„(t) 



(31) 



(32) 



(33) 



(34) 



«=o 



Finally, for the time variation of |0(0) at finite volume V, we 
obtain 



dt 



10(0) 



1 las 



2 \ dt 



10(0) 



(35) 



which shows that /(1/2) (dS/df) acts as t he g enerator of time- 
translations. As observed elsewhere it is remarkable 
that the same dynamical variable S whose expectation value is 
formally the entropy also controls time evolution: A privileged 
direction in time evolution (arrow of time) emerges which sig- 
nals the breaking of time-reversal invariance. 



4. Free energy and entanglement 

By acting with the operator J2 on the operators Uk and bk 
we obtain the Bogoliubov transformations ak 1-^ ak{0{t)), bk i-> 

bkiem 



akiOit)) = e''2 2qr,y,,,^^g-,Y2E,r,y,,, 

= dk cosh 0ii(t) - fo^ sinh 0^(0 , 

bk(0(t)) = e''2Zqr,y2,,^^g-,r2Z,r,y2,, 

= -a^ sinh 6kit) + bk cosh 6ii(t) 



(36) 



(37) 



where Okit) = F^f. The ak(6{t)) and bk(d{t)) operators are the 
annihilation operators for the vacuum |0(0) since 



flk(0(O)|O(O) = = M0(O)|O(O) , Vf. 

At each instant t and for each k, we have 

riaMt)) = <0(OI«I«k|0(0) = sinh2 0i(O . 



(38) 



(39) 



Similarly one obtains nt,^(6(t)) - sinh^ 0t(f) for the modes. 

We have observed that the state |0(f)) is an SU(1,1) general- 
ized coherent state. Eq. (39[ and the similar one for nh^(6(t)) 
show that it is a coherent condensate of equal number of at 
and bk modes. One can show that the creation of the ak(S(0) 
mode is equivalent to the destruction of the bk(0(t)) mode and 
vice-versa dUlt]. This means, the b\i,{0{t)) modes can be in- 
terpreted as the holes for the ak(^(0) modes and vice-versa. In 
other words, the /7(0(f))-system can be considered as the sink 
where the energy dissipated by the fl(0(f))-system flows and 
vice-versa. 

In this context, we also note that 



(40) 



Thus the difference Sa -Sb is constant under the time evolution. 
Since the fe(0(f))-modes are the holes for the fl(0(f))-modes, Sg- 
Sb is in fact the entropy for the closed system. 

By closely following Ref. |21] we introduce the free energy 
functional for the fl(0(f))-modes (we could do it as well for the 
^(0(f))-modes) 



Fa = {0(t)\(Ho,a - ^S,j0(t)) 



(41) 



where Ho^a = Zk^/t'^k'^k with Ek = Ak/2. Assuming that /3 - 
1 /ksTit) (ks denotes the Boltzmann constant) is a slowly vary- 
ing function of the time t, the stability condition dFgldOk - 0, 
Vit, gives pEk = - In tanh^ 0k{t). Then 



«„,(0(O) = sinh2 0i(f) 



1 



ePrnt - 1 



(42) 



which is the Rose distribution for ay at time t. We thus recog- 
nize that {|0(r))) is a representation of the CCR at finite temper- 
ature llisll. We can show that 



(43) 



dFa - dEk - ^dS = . 
Indeed we find 



dEk = J^Et hajm)dt = ^dS = dQ 



(44) 



Here hajfi(t)) denotes the time derivative of ria^iOit)), and, as 
usual, we define heat as dQ - {l//3)dS . 

We finally remark that the vacuum |0(f)) can be written in the 
following form 



io(f)> = n 



1 



cosh0jt(f) 



[|0)®|0> 



+ tanh 0t{t) (Iflk) ® l^k)) + ■■■]. 



This shows that it cannot be factorized into the product of two 
single-mode states: it is an entangled state for the modes Ok and 
b\i. Eqs. (I3ni-(l34li then show that S provides a measure of the 
degree of entanglement: the probability of having entanglement 
of the two sets of n a-modes and n b-modes is W„. Since W„ is 
a decreasing monotonic function of n, the entanglement is sup- 
pressed for large n. It appears then, that only a finite number of 



entangled terms in the expansion (ISTT i is relevant. However, this 
is only true at finite volume. The entanglement is truly realized 
in the infinite volume limit, i.e. in QFT, where the summation 
in Eq. OTb extends to an infinite number of components and 
Eq. dllll holds |27]. 

Notice that the robustness of the entanglement is rooted in the 
fact that, once the infinite volume limit is reached, there is no 
unitary generator able to disentangle the a and b modes. Such 
a non-unitarity is only realized when all the terms in the series 
( 31 i are summed up, which indeed happens in the V — > oo limit 
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5. Conclusions 

The discussion presented in this paper leads us to conclude 
that, provided the resonance condition co^ x does not ap- 
ply, the counter-rotating part Hqr of the Hamiltonian of two 
linearly coupled bosonic reservoirs shown in Fig. [T] may re- 
veal very interesting dynamical features which are typical of 
non-perturbative quantum field theory. By resorting to known 
results ilSl 12111 . the discussion of the explicit expression of 
(010(0) in such a frame for a concrete quantum optical system 
has shown that the RWA cannot be applied in conditions far ofi" 
resonance. There are absolutely non-trivial features in the phys- 
ical behavior of this system which have been overlooked so far 
to our knowledge. 

In particular, for non-vanishing Hcr turns out to be re- 
lated to the entropy operator This signals irreversible (non- 
unitary) time evolution, which is a manifestation of the break- 
down of time-reversal invariance (the arrow of time). Time evo- 
lution of the system shown in Fig.[T]hence needs to be described 
in terms of "trajectories" in the space of the unitarily inequiva- 
lent representations {|0(f))}- The free energy functional is mini- 
mized on such trajectories (at each time t we have dF - 0). In 
each representation, the system ground state |0(f)) turns out to 
be a generalized coherent state which is an entangled state of 
the modes in terms of which Hcr is expressed. 

Beyond their theoretical interest, our results may have some 
relevance also from the experimental standpoint 0151 For 
example, a study of a possible energy concentrating mechanism 
in atom-cavity system, which might become the object of an 
actual experimental observation, is currently in progress with 
preliminary results reported in Ref. jisll . 
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Appendix 

We present below some formulas used in the derivations in 
the text. For simplicity we omit the momentum suffices k. Let 



6 



us start by presenting the generators of the SU(1,1) group: 



Ji = ^(-/+ + J-) = ^(«'^' + ' 

Ji^-'-{J^-J-)^-'-{a'b'-ab), 

73 = ^(fl'fl + b^b + 1) , 
or else 

J+ - b^ , J- — ab , 
with commutators 

= -273, [73,^] = ±/±. 



(A.l) 
(A.2) 
(A.3) 

(A.4) 
(A.5) 



Notice that when the suffix k is introduced, we may define 7,- = 
Zk -^ki , ' = 1, 2, 3, and then the group structure is the one of 
SU{\,\)- S t/k(l , !)■ We also introduce the generators of 
the SU(2) group: 



h = ^(/+ + /-) = ^-{a'b + b^a) , 
h^-'-{h-I-)^-'-(a^b-b^a), 
h = ^(a^a - b'^b) , 

and 

1+ - a^b, I- - b'a , 
with commutators 

[/+,/-] =2/3, [/3,/J=+4. 



(A.6) 
(A.7) 
(A.8) 

(A.9) 
(A. 10) 



Notice that the Casimir operator for the 7's su(l,l) algebra is 73, 
while the Casimir operator for the /'s su(2) algebra is 73 - j = 
jia^a + b'^'b), i.e. 

[h-^Ji]^0, [h,Ji]^0, / = ±3. (A.ll) 

The K's and K's generators given below also close the su(l,l) 
algebra: 



Ki = ^[(a^ + a'^) + (b^ + b'^^)] , 
K2 = ^[(a^ - a'^) - (b^ - b"^^)] , 

K3 = ^[a'a-b"'b] ^ I3 . 
with commutators 

[KuK2]^-iKi, [^3,^1,2] = ±/^2,i 

and 

^1 = ^[(a^+a^^)-(b^ + b^^)], 
K2^^[(a^-a^'^) + (b^-b^'^)], 
K3 = ]^Wa-b^b] =Ki=h . 



(A. 12) 
(A.13) 
(A. 14) 

(A.15) 

(A. 16) 
(A. 17) 
(A.18) 



with commutators 

[KuK2] = -iK3, [K3,Ki,2] = +iK2,i . (A.19) 

K2 is the squeezing generator fl^. Note that the K's and K's 
can be constructed by using convenient combinations of the 
generators Kg-j , cr - a,b ; 1,2,3: Ko-,i - j(cr + cr' ), 

Ko-,2 - i(cr^-<T^\ Ka-,3 - 5(cr"'o-+ 5), which close the su(l,l) 
algebra for each cr. Other formulas used in the text are: 

2e (jj 
-4- , cos - 



sin0 : 



4^' 

sinha , cosh a : 

B 



A 
a>+A 



(A.20) 
(A.21) 



which have been obtained from Eq. (fT2l) and where 

A = yja? + 4^2 ^ B = yJojlA^ - 16^4 ^ (^.22) 

with B assumed to be different from zero: B 0. Moreover, we 
have also used 

[72, /i] = iKi , [72, K2]=0, [K2, h] = iKi , (A.23) 

[K2, h] = iKi , [K2, Ki] = iJ3 , (A.24) 

and obtained 



(A.25) 



w+(73-^)+ A 73-^^72 
-— ^1, 



H{e,a) ^ e'"'^m(e) e-'"'^^ 



(A.26) 



2gw_ 



B 1 

= —73 (jJ+ H h 

A 2 B A 
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